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Abstract. Here, we show that there is no positive integer n such that 
the nth Cullen number C„ = n2" + 1 has the property that it is com- 
posite but 4>{Cn) I C„ — 1. 

1. Introduction 

A Cullen number is a number of the form C„ = n2" + 1 for some n > 1. 
They attracted attention of researchers since it seems that it is hard to find 
primes of this form. Indeed, Hooley [8] showed that for most n the number 
Cn is composite. For more about testing C„ for primahty, see [3] and [6j. 
For an integer a > 1, a pseudoprime to base a is a compositive positive 
integer m such that a"^ = a (mod m). Pseudoprime Cullen numbers have 
also been studied. For example, in [12j it is shown that for most n, Cn is 
not a base o-pseudoprime. Some computer searchers up to several millions 
did not turn up any pseudo-prime C„ to any base. Thus, it would seem that 
Cullen numbers which are pseudoprimes are very scarce. A Carmichael 
number is a positive integer m which is a base a pseudoprime for any a. 
A composite integer m is called a Lehmer number if (j){m) | m — 1, where 
0(m) is the Euler function of m. Lehmer numbers are Carmichael numbers; 
hence, pseudoprimes in every base. No Lehmer number is known, although 
it is known that there are no Lehmer numbers in certain sequences, such as 
the Fibonacci sequence (see [9]), or the sequence of repunits in base g for 
any g G [2, 1000] (see j4]). For other results on Lehmer numbers, see [l], [2], 

m, m, m- 

Our result here is that there is no Cullen number with the Lehmer prop- 
erty. Hence, if (/>(C„) | C„ — 1, then C„ is prime. 

Theorem 1. Let Cn be the nth Cullen number. If <f){Cn) | — 1, then Cn 
is prime. 

2. Proof of Theorem [T] 

Assume that n > 30, that (j){Cn) \ Cn — ^, but that Cn is not prime. Then 
Cn is square-free. Write 



Cn = Y\.Pi 
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So, 

k 

- 1) I n2\ 

i=l 

Write n = 2"ni, where ni is odd. Then C„ = ni2"2 + 1, where 77-2 := a + n. 
Let p be any prime factor of Cn- Since p — 1 \ Cn — ^, it follows that 
p = mp2"p + 1 for some odd divisor nip of n and some Up with 

log n 

np<n2=n + a<n + -. 

log 2 

Let us first show that in fact Up < n. Assume that Up > n. Then, 

(1) C7„ = n2" + 1 = pA, 

for some positive integer A, where p > 2"+^ + 1. Observe that A > 1 because 
Cn is not prime. Now 

^ = V-2^TrTT<™- 

Reducing equation ([T|) modulo 2", we get that 2" | A — 1, so 2^^ < A — 1 < n, 
which is false for any n > 1. Hence, Up < n. 

Next we look at nip. limp = 1, thenp = 2*^^ + 1 is a Fermat prime. Hence, 
np = 2"^^ for some nonnegative integer 7. Since 2'^p = Up < n, we get that 
7p < (log n) I (log 2). Hence, the prime p can take at most 1 + (log n) /(log 2) 
values. Next, observe that since 

(2) JJ mp I n, 

p|C„ 

it follows that the number of prime factors p of Cn such that nip > 1 is 
< (log n)/ (log 3). Hence, we arrived at the bound 

/ s , logn logn 

3 A:<l + -^ + -^<l + 2.41ogn. 

log 2 log 3 



We next bound Up. Put := [y^n/log nj, and consider pairs (a, 6) of 
integers in {0, 1, . . . , A}. There are (A + 1)^ > nj log n such pairs. For each 
such pair, consider the expression L(a,6) := an + bup e [0,2n^/'^ /{logn)^/^]. 
Thus, there exist two pairs (a, 6) 7^ (ai,6i) such that 

2j^3/2 /(iQg 7^)1/2 

\(a—ai)n+(b—bi)np\ = \L(a,b)—L(ai,bi)\ < — < 3(n log n)"*^/^. 

n I log n — 1 

Put u := a — ai, v := b — bi. Then {u, v) ^ (0, 0) and 

\un + vnp\ < 3(n log n)^^'^. 

We may also assume that u and v are coprime, for if not, we replace the 
pair (n, by the pair {ui,vi), where d := gcd{u,v), Ui := u/d, vi := v/d, 
and the properties that max{|ni|, < {n/ logn)^^^ and \uin + viUpl < 
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3(nlogn)^/^ are still fulfilled. Finally, up to replacing the pair {u,v) by the 
pair (— n, —v), we may assume that u >0. 

Now consider the congruences n2" = —1 (mod p) and mp2"p = —1 
(mod p). Observe that 2, n, rup are all three coprime to p. Raise the 
first congruence to u and the second to v and multiply them to get 

n"m^2"^+"f^ = (-1)"+^ (mod p). 

Hence, p divides the numerator of the rational number 

(4) A := n"m;2""+"'''' - (-1)"+^ 

Let us show that A 0. Assume that A = 0. Recall that C„ = ni2"'^ + 1. 
Thus, expression (jH) is 

A = n^m;2("+")"+"''^ - (-1)"+'^ = 0. 

Then n^rrip = 1, (n + a)u + vup = 0, and u + v is even. Since ti > 0, it 
follows that V < 0. Put w := —v, so w > 0. There exists a positive integer 
p which is odd such that ni = p"^ and rUp = p". Since u and v are coprime 
and u + V is even, it follows that u and v are both odd. Hence, w is also 
odd. Also, since nip divides rii, it follows that u < w. We now get 

(2"^"" + a)u - WHp = 0, 

so 

u w 



Up 2"p"' + a 

The left-hand side is > u/n = u/(2"p"), because Up < n = 2°p". Hence, 
we get that 

u u w u w w 
< — = leadmg to — < -. — ; — - < — . 

2°p« ~ Up 2°'p'^ + a p« ~ p'" + (a/2°) ~ p^ 

For /9 > 3, the function s ^ s/p^ is decreasing for s > 0, so the above 
inequality together with the fact that u < w implies that u = w (so both 
are 1 because they are coprime), and that all the intermediary inequalities 
are also equalities. This means that u = w = l, a = and n = Up, but all 
this is possible only when C„ = p, which is not allowed. If p = 1, we then 
get that ni = 1, so every prime factor p of Cn is a Fermat prime. Hence, we 
get 

k 

C„ = 2"^ + 1 = 11(2^'"'' + 1) = Yl 2^""^ ' 

i=l /C{l,...,fc} 

and k > 2, but this is impossible by the unicity of the binary expansion of 

Cn- 

Thus, it is not possible for the expression A shown at ([4]) to be zero. 
The size of the numerator of A is at most 

2T'+\nu+npv\^u^\v\ ^ 2l+3("log»^)^''^^2(ra/logn)i/2 
^ 2l+3("log")^''^ + (2/log2)(nlogn.)i/2 ^ 26(n log n)!/^ 
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In the above chain of inequalities, we used the fact that 3 + 2/ log 2 < 5.9, 
together with the fact that (nlogn)^/^ > 10 for n > 30. Thus, for n > 30, 
we have that the inequality 

(5) ^<26("i°g")'^' 

holds for all prime factors p of C„. 
Thus, we get the inequality 

k k 
1=1 i=l 

leading to 

(6) k > 



n 



1/2 



6(log 

Comparing estimates dSj) and ([6]), we get 

„l/2 



< 1 + 2.4 log n, 



6(logn)^/2 

implying n < 6 x 10^. 

It remains to lower this bound. We first lower it to n < 93000. Indeed, 
first note that since n < 6 x 10^, it follows that ii p = = 2^^ + 1 is a 
Fermat prime dividing Cn, then 7 < 18. The only such Fermat primes are 
for 7 G {0,1,2,3,4}. Furthermore, (logn)/(log3) < log(6 x 10^)/(log3) = 
12.1104 . . . Hence, /c < 5 + 12 = 17. It then follows, by equation ([6]), that 

„l/2 

< 17, 



6(logn)i/2 

so n < 122000. But then (log n)/(log 3) < log(122000)/(log 3) = 10.6605 . . . , 
giving that in fact k < 15. Inequality ([6]) shows that 

^^^<15 
6(logn)i/2 ' 

so n < 93000. Next let us observe that if n is not a multiple of 3, then relation 
([2]) leads easily to the conclusion that the number of prime factors p of C„ 
with nip > 1 is in fact < (logn)/(log 5) = 7.15338 .... Hence, the number of 
such primes is < 7, giving that k < 12, which contradicts a result of Cohen 
and Hagis [5] who showed that every number with the Lehmer property must 
have at least 14 distinct prime factors. Hence, 3 | n, which shows that C„ is 
not a multiple of 3. An argument similar to one used before proves that n 
is not a multiple of any prime q > 3. Indeed, for if it were, then relation ([2]) 
would lead to the conclusion that the number of prime factors p of C„ with 
mp > 1 is < 1 + log(n/g)/(log3) < 1 + log(93000/5)/(log 3) = 9.94849..., 
so there are at most 9 such primes. Also, C„ can be divisible with at most 
4 of the 5 Fermat primes with 7 G {0,1,2,3,4}, because 3 = Fq does 
not divide Cn- Hence, /c < 9 + 4 = 13, which again contradicts the result 
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from [5j. Thus, n = 2" 3^ and so all prime factors p of C„ are of the form 
2"i3/3i _)_ 1 for some nonnegative integers ai and Now write 



for some integer a > 2. Since 




n (^"^2°i3/'i) 

0, /3i>0 ^ ^ 
1+1 prime 



we get that a < 2, which is a contradiction. This shows that in fact there 
are no numbers Cn with the claimed property. 

We end with some challenges for the reader. 

Reearch problem. Prove that Cn is not a Camichael number for any 
n > 1. 

If this is too hard, can one at least give a sharp upper bound on the count- 
ing function of the set C of positive integers n such that C„ is a Carmichael 
number? We recall that Heppner [7] proved that if x is large then the num- 
ber of positive integers n < x such that Cn is prime is 0{x/ \ogx), whereas in 
[12j it was shown that if o > 1 is a fixed integer then the number of positive 
integers n <x such that Cn is base a-pseudoprime is 0(x(log log x)/ log x). 
Clearly, imposing that Cn is Carmichael (which is a stronger condition) 
should lead to sharper upper bounds for the counting function of such in- 
dices re. 

Finally, here is a problem suggested to us by the referee. Theorem [T] shows 
that 4>{Cn)/ gcd{Cn — 1, 4>{Cn)) exceeds 1 for all re. Can one say something 
more about this ratio? For example, it is possible that a minor modification 
of the arguments in the paper would show that this function tends to infinity 
with re, but we have not worked out the details of such a deduction. It would 
be interesting to find a good (large) lower bound on this quantity which is 
valid for all n and which tends to infinity with n. How about for most n? 
What about lower and upper bounds on the average value of this function 
when re ranges in the interval [l,x] and x is a large real number? We leave 
these questions for further research. 
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